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Abstract 

By expanding the projection operator in powers of the density fluc- 
tuations, we conjecture a hamiltonian purely quadratic in the sea-bosons 
that reproduces the right spin and charge velocities and exponent for the 
U — oo case in one dimension known from the work of Schulz. Then 
we argue that by simply promoting wavenumbers to wave vectors we are 
able to study the two dimensional case. We find that the quasiparticle 
residue takes a value Zf = 0.79 close to half-filling where it is the small- 
est. This is in exact agreement with the prediction by Castro-Neto and 
Fradkin nearly ten years ago. We also compute the magnetic suceptibility 
and find that it diverges close to half-filling consistent with Nagakoka's 
theorem. 

1 Introduction 

The large-U Hubbard Model or the t- J model and its variants have been the sub- 
ject of active study ever since its importance to understanding cuprates has been 
realised. For a review see the article by Dagotto[l]. The slave-boson/fermion 
approach that naturally takes into account the feature of spin-charge separa- 
tion has been employed by various authors [2] (for some recent references we 
note the work of Balents et.al.[3] and Wang[4]). Recently the phase diagram of 
the t-J-V model has been found by using a thorough but tedious procedure of 
linked cluster expansion by Zheng et.al.[5]. High temperature expansions have 
been used to study the t-J-V models and the t-J models by various groups [6] 
[7]. The cluster dynamical mean-field approach has been employed by Jar- 
rell et.al.[8] to compute the phase diagram of the Hubbard model. We wish 
to complement their study and others like them through a simpler analytical 
approach. The case of the infinite U Hubbard model has been considered by 
Nagaoka who shows that the ground state is a ferromagnet [9] for a single hole in 
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an otherwise half-filled band. Shastry, Krishnamurthy and Anderson [10] have 
shown that this ground state is unstable for large enough hole concentration. 
Gutzwiller projected variational wavcfunctions[ll] have also been used to study 
the t-J model. Chen and Tremblay[12] have used Monte Carlo simulations to 
compute the magnetic susceptibility for the large but finite U version. Mishra 
and Kishore and Mishra[13] have studied the thermodynamics of the infinite 
U Hubbard model and the issue of spin-charge separation using the method of 
orthofermions. We shall address the issue of spin-charge separation when the 
author's hydrodynamic approach[14] to the Hubbard model is accepted by the 
community. 

In this article, we wish to write down an effective low energy theory of the 
t-J=0 model that is local in the operators in the sea-boson language (RP A- type) 
that involves introducing 'renormalised' doping dependent hopping and onsite 
terms. Thus the claim is that one may treat the infinite U Hubbard model the 
same way as we treated the small U Hubbard modcl[19](RPA-like) provided we 
'renormalise' the parameters in the hamiltonian. We find that this leads to non- 
trivial predictions for the quasiparticle residue and velocity of the quasiparticles 
in more than one dimension. We also compute the magnetic suscpctibility and 
show that it diverges near half-filling indicating that the system exhibits fer- 
romagnetic instability arbitrarily close to half-filling, consistent with Nagaoka's 
theorem. This also indicates that the formula for the residue is also reliable as 
is the conjuctured low energy hamiltonian for the t-J model that is local in the 
operators. This then sets the stage to study the large but finite U Hubbard 
model and also the t-J-V models in future publications . 



2 The Theory 

Here we describe the theory which we are going to use. Consider the t-J model 
in one dimension with U = oo. 




r = H(l-n^ nil ) (2) 

i 

We would now like to recast this in the sea-boson language so that we may 
recover the results of Schulz[15]. For this we mentally expand the projection 
operator in powers of the density fluctuations and retain the leading terms. 
Thus we could write for example, 

■p _ e ^2 i Log(l-<n iT nn>+\ [nit ™»J.-<™tt ™U>]) ^ 

and expand in powers of A = — 1 and retain the leading terms. We shall not 
do this explicitly however we may redefine effective parameters that simulate 
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such an expansion. The idea is to arrive at useful answers quickly and effort- 
lessly. The above procedure for example, entails the solution of self-consistent 
equations for the density correlation functions, this is clearly not desirable as 
it is likely to be complicated. We adopt the point of view that making contact 
with the Id system allows us to generalise to the case of more than one dimen- 
sion by simply promoting wavenumbers to wavevectors. This we justify mainly 
by pointing out that the magnetic susceptibility derived using this approach is 
consistent with the rigorous Nagaoka theorem[9]. 

This means we may suspect that the t- J model with J = also has the form, 

H = e k c\ a c k<J + ^jti- PqlP-ql (4) 

for some suitable 'renormalized' dispersion e k = —Ueff cos(ka) and U e ff. The 

i 

spin and charge velocities may be written as follows. v s — vp, e ff — ^J S SS ) ^ 

and v c = vp,eff + ^J s fS ) 2 > where VF, e ff — ^eff sin(ir n e /2) (lattice spac- 
ing a = 1). From the work of Schulz we know, v c = 2t sin(ir n e ). Here 
n e — N e /N a is the number of electrons per site and v s — 0. Also we know 
that the anomalous exponent for the momentum distribution is 7 = 1/8. Thus 
we would like to compute t e // and U e ff in terms of t and n e given these 
facts. From our earlier work[19] we may read off a formula for the exponent 

: 7 - U & ,V %" K+^TTF " Thus we ma y write > u eff - t ^/^TJ ("c + 
VF.eff) and vp e ff + VF.eff — w c — we m ay solve this to yield, VF.eff = 
+ \J + 4«c^ • We choose the positive solution, U e ff > since we 

know that the negative U Hubbard model possesses a spin gap [17] whereas the 
U = 00 Hubbard model docs not. These may be solved via a scaling argument, 
U e ff — u e ff v c and VF, e ff = Ueff v c - This may be solved using mathemat- 
ica to yield u e ff = 4.71 and y e ff = 0.5. In order for v s = we must have 
U e ff > irVF, e ff- We may see that this is being obeyed. Hence we may claim 
that Eq.( 4) captures the low energy physics of the Hubbard model at U = 00 
exactly. Now we make the following non-obvious assertion that the same Eq. ( 4) 
also captures the same physics in more than one dimension by simply promoting 
wavenumbers to wavevectors. We make make this plausible by pointing out that 
the one dimensional Hubbard model may be generalised to higher dimensions by 
precisely such a procedure, namely by promoting wavenumbers to wavevectors. 
On the other hand, the prediction (see below) that the magnetic sucseptibility 
diverges as half-filling is approached is consistent with the rigorous Nagaoka 
theorem, thus lending credibility to our approach. 
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2.1 Zp in Two Dimensions 



The t-J model in two dimensions is relevant to cuprates[2] [11] and cobalt oxide 
superconductors[16]. In the two-dimensional case, we may solve for the sea- 
boson occupation as follows [19]. The energy dispersion for a square lattice is 
£k = —2t e ff [cos(k x a) + cos(k y a)]. The sea-boson method is well-known and 
hence we shall refer the reader to our earlier work for details[19]. Suffice it to 
say that we have to compute the boson occupation numbers in order to derive 
a formula for the quasiparticlc residue. Therefore we write, 



< 



N a e c( C li ^c) {-U) c - 6 k + £k-qj 



(5) 



< 4^(* + )4£ q * 



u, 



eff 



[A?_ 



q5-» ^k-qffJ 



N a e,(q, u s ) (-w s - e k + £k- q ) 2 



(6) 



where P(q, W ) = £ k ^^fg^) and £c(qjC ,) = i _ £|/ P 2 (qjW)j alsQ; 

e s (q, w) = 1 + % £i P(q, w). In order to evaluate these experssions, we turn 
the cosine dispersion into a parabolic one by demanding that the slope of the 

,2 

two be the same at the Fermi momentum. —2t e ff cos{k x ) — —2t e ff + 
Thus we set 2t e ff sin{kp /^/2) = kp/(\/2m). Also, kp = \J2-k n e . This approx- 
imation captures the important physics while leaving the integrals analytically 
computable. From our earlier work we may read off the formula for the velocity 
of the charge carriers. Here vp,2d — kp/m. We find that there is only one 
velocity and e s = does not have a solution. 



1 + 



VF.2d 



(7) 



and the formula for the quasiparticle residue has been derived in an earlier 
work[19]. 

Close to half-filling, 2t e jf w y e fj 2t sin(ir n e ), U e ff — u e ff 2t sin{-K n e ). 
Furthermore, y e ff 2t sin(n n e ) sin(y/Tr) = v^m- Therefore we may write the 
following formula for the quasiparticle residue. 



Zp = Exp 



2 \4 



irkpinvc 



( 



Amv cc ArcTan 



TT — 



mvcc — k£ 



m 2 v' 2 < , —k 2 F 



k'p 



(8) 



This may be evaluated to yield, 



0.79 



(9) 
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This result is in exact agreement with the prediction by Castro-Neto and Fradkin 
nearly ten years ago [18] . They considered spinless fermions interacting via short- 
range interactions in two dimensions. Thus it would appear that the same 
physics is operating. In general we may conclude that the same mechanism 
that makes the anomalous exponent saturate in one dimension to a rather small 
value (1/8) also makes the quasiparticle residue to saturate to an equally small 
deviation from unity. In one dimension it suggests that the residual Fermi 
surface persists all the way up to infinite onsite repulsion. In two dimensions 
it says that the actual Fermi surface persists all the way upto infinite onsite 
repulsion. Since the functional dependence on the onsite repulsion is monotonic, 
we may conclude that there is no chance for Fermi liquid theory to break down 
in two dimensions with short-range interactions. In fact if anything these results 
show that the system is strongly metallic all the way upto infinite repulsion. 

3 Magnetic Susceptibility 

The dynamic spin susceptibility may be computed using the Kubo formula. 

/•oo 

xH=i/ dte^ 1 <[5+(t),5-(0)]> (10) 
Jo 

s + W = E c k T (*Hi(t) (ii) 

k 

S~(0) = E c ki(0)ck T (0) (12) 

k 

To use the sea-boson formalism, we have to take special care to account for 
possible infrared divergences. This is the crucial aspect that leads to the correct 
solution of the Luttinger model[20] and also leads to the solution of the problem 
of quenched disorder for scattering across the Fermi surface[?]. Define nn(k) = 
c kT Ck l' Now we would like to compute the correlation function, 

JV(kt;kY) = (n n (k,i)n iT (k',i')} (13) 
This may be decomposed as follows. 

N(kt; kY) = (1 - n F (k))(l - n F (k '))S AA {kt; kY) + n F (k)n F (k')S , Bs(kt; kY) 
-(1 - n F (k))n F (k')S AB (kt;k't) - (1 - n F (k ))n F (k)S BA (kt; kY) (14) 

S mn (kt;kY) =e -2<S m , n( k)> e -2<s„, iT (k')> S° mn (kt;kt) (15) 
S A {kaa',t) = ^4_ q/2(7i (qa,i)A k _ q/2CTl (qY,i) (16) 
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S B (kaa , t) = ^ A [ +<l/2a > fa^i . ^k+q/ 2( x (q^i , *) 



(17) 



SS, n (k*;kY) =< S rn (kTi,t)S„(k' IT,*') > - < Sm(kU,t) >< S„(k' IT,*') > 

(18) 

Since the (a, <r) part of the hamiltonian is distinct from the (a, a) part of the 
hamiltonian, we may conclude that < Sa >~< Sb >= 0. Thus in this case 
there are no difficulties associated with the infra-red regulator. After some 
tedious calculations we may see that the most divergent part of the magnetic 
susceptibility is given by 

x{uj = 0) - m ~ (1 - n e ) _1 (19) 

As pointed out before, this is consistent with Nagaoka's theorem and also with 
Shastry et.al.'s[10] work that shows that the ferromagnetic ground state is un- 
stable with respect to the addition of holes. We find that unless the filling is 
arbitrarily close to half-filling, the susceptibility does not diverge. 

4 Conclusions 

To conclude, we have computed the quasiparticle residue of the infinite-U Hub- 
bard model in two dimensions saturates to a value Zf = 0.79. This is in exact 
agreement with the prediction of Castro-Neto and Fradkin[18]. The velocity of 
the quasiparticles is also found to shrink to zero as expected. The magnetic 
susceptibility is found to diverge as half-filling is approached which is consistent 
with a ferromagnetic ground state for a small 1 concentration of holes. This 
is consistent with Nagaoka's theorem[9] and also the with the work of Shastry 
et.al.[10]. This shows that the effective low energy hamiltonian for the t-J model 
that is local in the operators is now reliable and can be expected to yield similar 
nontrivial results when generalised to include J terms and so on. 

It is a pleasure to acknowledge email correspondence with Debanand Sa and 
Pinaki Majumdar and useful comments by Sourin Das. 
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